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ABSTRACT 


Rayleigh  wave  group  velocities  were  obtained  for  some 
Eurasian  paths  using  VLPE  data  for  events  occurring  during  the  ISM. 

Fundamental  Rayleigh  wave  group  velocities  and  spectra  were 
determined  using  narrowband  filters  to  minimize  possible  multipath  and 
higher  mode  effects.  Further,  the  spectra  were  corrected  tor  geometrical 
spreading  and  effective  attenuation. 

Observed  Rayleigh  wav  group  velocities  for  propagating  paths 

c *ossing  Central  India,  and  the  Russian  and  Siberian  platforms  were  in 
close  agreement  to  those  observed  for  the  Central  United  States.  The 
observed  group  velocities  and  amplitude  spectra  for  paths  in  and  normal  to 
the  Alpine -Himalayan  Fold  System  indicate  the  presence  of  thick  crustal 
sect  ons  having  high  effective  attenuation  properties. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and  this 
document  is  subject  to  later  revision  as  may  be  necessary.  The  views  and  con- 
clusions presented  are  those  of  the  authors  and  should  not  be  interpreted  as 
necessarily  representing  the  official  policies,  either  expressed  or  implied,  ot 
the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  US  Government. 
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SECTION  1 
INTRODUCTION 

The  establishment  of  the  Very  Long  Period  Experiment  (VLPE) 
provided  a small  netwo-k  of  high  gain,  high  quality,  long-period  digital  seismo- 
graphs at  var'ous  locations  throughout  the  world.  Further,  the  availability 
of  a list  of  confirmed  and  well  located  events  occurring  during  the  International 
Seismological  Month  (ISM)  provided  the  opportunity  to  determine  and  compare 
Rayleigh  wave  group  velocities  of  several  continental  travel  paths  in  Eurasia  to 

those  determined  by  others. 

For  this  study,  we  selected  three  central  Asian  events  which 
occurred  during  the  ISM  and  were  recorded  at  several  VLPE  stations  located 
in  Eurasia.  Fundament?  I Rayleigh  wave  group  velocities  and  spectra  were 
determined  using  narrowband  filters  to  minimize  possible  noise,  multipath, 
and  high  mode  effects. 

The  spectra  were  corrected  for  geometrical  spreading  and 
effective  attenuation.  Turnbull,  et  al.  (1973)  using  Tryggvason's  method  (1965) 
determined  the  effective  attenuation  for  a selected  path  from  Sinkiang  to  Chang 
Mai,  Thailand  (CHG). 

The  group  velocities  are  compared  to  those  determined  by 
Santo  and  ^at6  (1^66)  for  various  tectonic  regions  in  Eurasia  and  to  those 
determined  by  rtrunc  and  Dorman  (1963)  and  McEvilly  0q64)  for  the  mid- 
continental United  States. 

These  data  provide  the  basis  of  the  discussion  or  our  spectral 
estimates  and  group  velocities. 
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FFCTTON  II 

DATA  AND  METHOD  OF  ANALYSIS 


A.  DATA  USED 

Seismic  recordings  from  the  VLPE  stations  at  Chiang  Mai, 
Thailand  (CHG),  Fairbanks,  Alaska  (FBK),  Toledo,  Spain  (TLO)  and  Kongsberg, 
Norway  (KCN)  were  selected  for  three  earthquakes  located  near  the  I SSR- 
Mongolia  border,  in  Tadzhik,  and  southwestern  Afghanistan.  Event  and 
station  parameters  are  given  in  Table  I.  From  these  twelve  station-event 
combinations  we  have  long-period  data  for  ten  of  the  possible  Eurasian  paths. 

No  record  was  available  at  TLO  for  event  231*and  the  signal  to  noise  ratio  was 
too  small  for  any  type  of  signal  ana'ysis  of  the  record  at  FBK  for  event  392. 
Station  and  event  locations  along  with  the  great  circle  path  representations  are 
shown  in  Figure  II-L 


» V 


B.  METHOD  OF  ANALYSIS 

The  determination  of  fundamental  mode  Rayleigh  wave  group 
velocities  and  amplitude  spectra  was  accomplished  with  the  use  of  a series  of 
narrow  band  filters.  Briefly,  for  this  analysis  we  filter  in  the  frequency 
domain  with  a series  of  narrow  band  filters  with  specified  center  frequencies. 

The  center  frequencies  selected  correspond  to  periods  of  15  through  60  seconds 
at  5 second  increments  (i.  e.  T = 15,  20,  25,  . . .etc.).  The  inverse  Fourier 
transforms  of  these  computations  result  in  series  of  time  traces  having  well 
behaved  oscillating  functions.  The  specifics  of  the  method  are  given  by 
Alexander  (1963),  Archambeau  et  al  (1965),  and  recently  by  Turnbull  et  al  (1973). 
The  resultant  oscillating  wave  packets  correspond  to  group  energy  packets  with 
the  largest  peak  amplitude  generally  corresponding  both  in  time  and  amplitude  to 
the  dominant  mode.  Figure  II-2  is  a typical  example  of  the  application  of 
narrow  band  filtering. 
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TABLE  II-l 

EVENT  AND  STATION  PARAMETERS 
EVENT  PARAMETERS 


Event 

Number 

Date 

231 

2/26/72 

375 

3/17/72 

392 

3/20/72 

Area 


Lat. 

Long. 

O.  T. 

mb 

50.  6N 

97.  1 E 

23:31:04 

5.2 

40.  IN 

70.  0 E 

09:17:11 

5.2 

30.  ON 

60.  0 E 

20:08:31 

3.4 

VLPE  STATION  PARAMETERS  AND  EP1CENTRAL  DISTANCES 

| VLPE  Station  Parameters  Distance  in  Degrees 

Station  Designator  Lat. Long.  231  375  392 — 

Chiang  Mai.  Thailand  CHG  18.79N  98.  98  E 30.4  32.1  36.9 

Fairbanks,  Alaska  FBK  64.  90  N 148.  01 W 55.5  72.5  81.8 

Toledo,  Spain  TLO  39.  86  N 4.02V/  - ^4.  7 52.6 

Kongsberg,  Norway  KON  59. 65  N 9.  59  E 48.3  i 42.3  44.9 


STATION  AND  EVENT  LOCATIONS 


We  believe  the  determination  of  group  velocities  and  spectra 
by  this  method  to  be  intrinsically  more  reliable  than  other  methods  for  the 
following  reasons: 

• It  is  well  known  that  lateral  variations  along  the  propagating 
path  car  cause  deviations  from  the  expected  direction  of  the 
various  group  arrivals  of  surface  waves  (Evernden,  1953, 
Alexander,  1963,  Glover  and  Alexander,  1969,  and  Newton, 1973). 
These  effects  may  be  termed  as  multipathing. 

• Fourier  spectral  estimations  of  the  total  signal  time  series  may 
not  be  valid  for  the  fundamental  mod  since  some  earthquakes 
do  excite  significant  higher  mode  energy. 

• Bursts  of  noise  in  the  signal  gate  can  be  better  discerned  and 
discarded. 

Thus,  some  of  the  effects  of  noise,  multipathing  and  the  presence  of  higher 
modes  can  be  mh.’mized  with  the  proper  use  of  narrow  band  filtering  in  the 
determination  of  fundamental  Rayleigh  wave  group  velocities  and  spectra. 
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SECTION  III 
RESULTS 


A.  GROUP  VELOCITIES 

Fundamental  Rayle.gh  wave  group  velocities  for  ten  Eurasian 
paths  are  shown  in  Figure  II1-L  In  the  determination  of  these  group  veloci- 
ties, we  did  not  correct  for  either  the  source  phase  or  instrumental  phase 
response.  However,  the  epicentral  distances  are  large  and  the  V LPF 
instrument  amplitude  response  is  nearly  flat  between  periods  of  25  to  45 
seconds.  Thus,  the  errors  due  to  either  source  or  instrument  phase  terms 

should  be  insignificant. 

The  observed  g-oup  velocities  are  indicated  by  the  various 
defined  symbols  on  the  graphs.  Included  for  comparison  (solid  line)  is  the 
group  velocity  curve  determined  from  McEvilly's  (1964)  model  for  the  mid- 
continental United  States  crust  and  upper  mantle  structure.  This  model 
varies  only  slightly  from  Brune  and  Dormans  (1)63)  Canadian  Shield  (CANSD) 

model  (lable  III- 1 ) . 


Events  231  and  375  lie  on  the  northern  portions  of  massive 
fold  systems  and  392  is  located  on  the  southern  side  of  the  Alpine-Hima’ayan 
fold  system.  Paths  l through  4 from  events  231  and  375  to  FBK  and  KON, 
Figure  II- 1,  cross  the  Siberian,  Northeastern  Siberian,  .Vestern  Siberian 
and  Russian  platforms.  The  general  tectonic  regions  have  been  defined  by 
Rodriguez  (1969).  Further,  paths  5 and  6,  from  event  392  to  TLC  and  CHG 
cross  Iran,  Greece,  Turkey,  and  the  Mediterranean  region,  and  the  northern 
part  of  India  all  of  which  are  essentially  platform  regions.  The  observed 
group  velocity  curves  for  the  first  six  paths  crossing  platform  regions 
compare  closely  to  that  of  McEvdly's  mid- continental  United  States  and 


III- 1 


TABLE  III-l 

COMPARISON  OF  CANSD  AND  MCEVILLY'S  EARTH 
MODELS  FOR  THE  CENTRAL  UNITED  STATES 


. CANSD  (1963) 

MCEVILLY  (1964) 

H 

a 

fi 

P 

H 

a 

fi 

P 

6.  0 

5.  64 

3.  47 

2.  70 

11.  0 

6.  1 

3.  50 

2.  70 

10.  5 

6.  15 

3.  64 

2.  80 

9.  ) 

6.  4 

3.  68 

2.  90 

18.  7 

6.  60 

3.  85 

2.  85 

18.  0 

6.  7 

3.  94 

2.  90 

80.  0 

8.  10 

4.  72 

3.  30 

24.  0 

8.  15 

4.  75 

3.  30 

100.  0 

8.  2 

4.  54 

3.  44 

40.  0 

8.  2 

4.  61 

3.  30 

100.  0 

8.  3 

4.  51 

3.  53 

18.  0 

8.  2 

4.  45 

3.  30 

80.  0 

8.  7 

4.  76 

3.  60 

- 

8.  7 

4.  80 

3.  60 

H = Thickness  (km),  a - Compres sional  Wave  Velocity, 
/3  = Shear  Wave  Velocity,  P - Density 
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CANSD  dispersion  curve.  Other  observational  evidence  by  Santo  and  Sato 
(1966)  verify  our  group  velocity  determinations  for  platform  regions  in 
Eurasia. 

Paths  7 and  9,  Event  375  to  TLO  and  Event  392  to  KON,  are 
considerably  more  complex  in  that  they  include  the  Caspain  and  Black  Sea, 
crustal  structure  and  portions  of  the  Alpine-Himalayan  fold  system.  For 
these  paths  we  observe  lower  average  crustal  velocities  than  for  the  platform 
regions  (see  Figure  III-l). 

Path  8 (Event  375  to  CHG)  is  parallel  to  and  in  the  Alpine- 
Himalayan  fold  system  and  we  observe  an  overall  lower  group  velocity  < arve 
wh.ch  is  indicative  of  a substantial  increase  in  crustal  thickness  over  the 
travel  path  as  compared  to  the  platform  regions. 

Path  10,  event  231  to  CHG,  crosses  the  Central  Asian,  Central 
China  fold  systems  and  a portion  of  the  Tibet  platform.  We  were  unable  to 
determine  valid  group  velocities  for  this  path  and  from  inspection  of  the 
signal  time  series  and  narrow  band  output  in  Figure  III-2  it  is  clearly  evident 
that  multipathing  does  occur.  The  group  velocities  shown  in  Figure  IIJ-1 
were  determined  by  averaging  the  time  of  one  or  more  of  the  wave 
packets  from  the  narrow  band  filtering.  The  results  indicate  low  group 
velocities  over  a large  period  band  and  imply  large  crustal  thicknesses  of 
relatively  constant  velocity.  The  multipathing  is  more  likely  due  to  th  * 
extremely  complicated  path  which  is  comprised  of  the  two  major  fold  systems 
and  the  Tibet  platform. 

B.  GROUP  SPECTRA 

Figures  III-3,4and  5 show  the  group  spectra  by  event  (231, 

375,  and  392), corrected  for  instrument  response,  geometrical  spreading, 
and  effective  Q.  AH  spectra  were  corrected  according  to  Tryggvasou' s 
(1965)  effective  energy  dissipation  (k^)  curves  shown  in  Figure  III- 6. 
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Event  392 


I I 1 

20  30  40 

Period  (Seconds) 


FIGURE  III— S 

GROUP  SOURCE  SPECTRA  FOR  EVENT  392 


Turnbull,  1973 


ENERGY  ATTENUATION  COEFFICIENT,  k , FOR  EURASIAN  PATHS 


This  curve  was  derived  from  WWSSN  station  records  of  four  presumed 
Novaya  Zemlya  explosions  and  several  natural  earthquakes.  Turnbull  (1973) 
also  determined  in  a similar  manner  the  effective  energy  dissipation  from 
VLPE  records  of  two  presumed  Eastern  Kazakh  explosions.  His  results  are 
shown  with  the  open  circles  and  are  in  close  agreement  with  Tryggvason's 
results.  Further,  Turnbull  calculated  the  effective  energy  dissipation  for 
the  path  from  Eastern  Kazakh  to  CHG  (solid  dots).  This  path  crosses 
the  Alpine-Himalayan  Fold  System.  The  data  in  this  figure  provides  the  basis 
of  the  discussion  for  our  spectral  estimations. 

We  can  make  the  following  general  comments  concerning  the  source 
amplitude  spectra  shown  in  Figure  3,  4,  and  5. 

• The  spectral  amplitudes  for  all  events  with  the  exception  of 
the  look  direction  represented  by  CHG,  are  overlapping  and 
relatively  close  in  value  to  each  other.  These  differences  can 
be  attributed  to  the  event  radiation  patterns  and  to  the  lack 

of  precise  effective  energy  dissipation  'kp)  factors  for  each 
path. 

• The  souice  spectra  for  the  look  direction  represented  by  CHG, 
as  previously  stated,  was  corrected  using  Tryggvason's  k^. 

These  amplitudes  are  significantly  lower  than  those  determined 
from  the  other  stations,  e specially  at  the  longer  periods. 
Correcting  these  spectra  with  the  k^  determined  by  r"urnbull 
(Figure  6)will  result  in  increasing  the  amplitudes  by  a factor 
of  about  2.  The  shapes  of  the  spectra  will  not  change  apprec- 
iably. The  amplitude  levels  for  path  6 now  become  comparable 
to  the  others.  However,  for  paths  8 and  10,  only  the  amplitudes 
at  short  periods  (15  to  30  second  periods)  become  comparable 
while  for  the  longer  periods  the  amplitudes  remain  low. 
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SECTION  IV 

DISCUSSION  OF  RESULTS 

The  group  velocities  measured  over  the  range  of  periods  from 
15  to  60  seconds  are  sampling  the  crust  and  upper  mantle  structure  to  depths 
up  to  300  kilometers.  Thus,  differences  in  the  observed  group  velocity  curves reW 
ativeto  CANSDor  McEvilly's  reference  group  velocity  curve  can  be  explained  by 
variations  in  thicknesses  and  elastic  parameters  of  the  crust  and  upper  mantle. 
Figure  1 shows  typical  cross  sections  of  the  crust  and  upper  mantle  for  the 
appropriate  paths  used  in  this  study  (Rodriguez,  1969). 

The  observed  group  velocities  for  paths  1 through  6 (paths 
through  Platform  regions)  show  slightly  lower  velocities  at  periods  of  50  to 
60  seconds  than  McEvilly.  This  suggests  that  the  low  velocity  zone  in  the 
upper  mantle  is  thicker  than  the  indicated  58  kilometers  by  McEvilly.  It  is 
probable  that  the  thickness  of  200  kilometers  shown  for  the  CANSD  model 
would  have  provided  more  appropriate  velocity  values  for  these  longer  periods. 

The  observed  group  velocities  for  path  7,  which  crossed  the 
Caspian  and  Black  Sea  regions  reflects  the  thickening  of  the  low  velocity 
surfacial  layer,  thinning  of  the  crust  and  lower  moho  velocity.  Here,  the 
short  period  velocities  are  significantly  lower,  the  slope  for  the  intermediate 
periods  is  comparable,  and  the  overall  velocity  curve  is  lower  than  that 
observed  for  the  platform  regions.  The  cross  section  in  Figure  IV-1  for  the 
Caspian  and  Black  Sea  clearly  supports  these  general  observations. 

Path  9 is  a complicated  path  which  crosses  the  Alpine- 
Himalayan  fold  system,  Caspian  Sea  region  and  the  Russian  platform.  We 
observed  significantly  overall  lower  group  velocities  but  little  or  no  unusual 
variations  in  amplitude  levels.  Since  the  observed  group  velocities  are  aver- 
ages of  the  structures  traversed,  we  can  only  suggest  that  the  Alpine- 


Himalayan  fold  system  portion  of  the  path  was  small  in  distance  compared  to 
the  rest  of  me  path.  This  caused  lower  group  velocities  but  had  little  effect 
on  the  effective  attenuation  of  the  Rayleigh  wave. 

Paths  8 and  10  are  extremely  complex  in  that  path  8 traverses 
the  Alpine -Himalayan  fold  system  and  path  10  crosses  the  Central  Asian. 

Central  China  fold  systems  and  a portion  of  the  Tibet  platform.  The  fold 
systems  result  f 'om  the  collision  of  continents  according  to  the  theory  of 
plate  tectonics  (Metz  and  Hammond,  1974).  The  Himalayas  are  thought  to 
have  termed  when  the  Indian  subcontinent  collided  with  the  Eurasian  contment. 
Little  is  known  about  mountain  building  caused  by  continental  collisions  since 
vertical  motions  are  almost  as  rapid  is  horizontal  motions. 

Thus,  the  cross  section  shown  for  the  Alpine-Himal  yan  and 
for  the  other  fold  systems  (Figure  IV-1)  is  an  over  simplification  of  the  true 
structure.  However,  it  does  provide  an  explanation  for  the  observed  low 
group  velocities  and  amplitudes  by  indicating  large  crustal  thicknesses  totaling 

about  75  kilometers. 
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FIGURE  IV-1 


TYPICAL  CROSS  SECTIONS  FOR  PLATFORM,  CASPIAN 

and  black  sea,  and  fold  systems 


SECTION  V 
SUMMARY 


Eurasia  is  comprised  of  vast  platforn  structures.  The  observed 
Rayleigh  wave  group  velocities  indicate  that  the  structure  of  the  crust  and  upper 
mantle  for  these  regions  is  similar  to  that  of  the  Central  l nited  Mates, 
ous  authors  have  written  that  all  continents  have  essentially  the  same  crust 
and  upper  mantle  structures.  Thus,  these  results  were  as  expected. 

However,  for  paths  encompassing  fold  systems  such  as  the 
Alpine -Himalayan,  Central  China,  and  Central  Asian,  we  found  a significant 
decrease  in  the  Rayleigh  wave  group  velocities  and  amplitudes,  indicating  high 
effective  attenuation  properties.  Further,  for  path  10,  multipathi  g is  evident. 


V-l 


SECTION  VI 
ACKNOWLEDGMENTS 

The  author  wishes  to  thank  T.  W.  Harley  and  L. 
for  their  helpful  discussions  and  critical  reading  of  this  paper, 
grateful  to  J.  S.  Shaub  for  his  capable  assistance  in  processing 
of  the  data. 


S.  Turnbull 
I am  also 
and  reduction 


Vl-l 


Alexander,  S.  S.  , 1963,  Surface  Wave  Propagation  in  the  Western  United 

States,  Ph.  D.  Thesis,  Calif.  Inst,  of  Tech.  , Pasadena,  California. 

Archambeau,  C.  B.  , J.  C.  Bradford,  P.  W.  Broom,  W.  C.  Dean,  E.  A. 
Flinn,  and  R.  L.  Sax,  1965,  Data  Processing  Techniques  for  the 
Detection  and  Interpretation  of  Teleseismic  Signals,  Proc.  IEEE, 

53,  1860  - 1884. 

Brune,  J.  and  J.  Dorman,  1963,  Seismic  Waves  and  Earth  Structure  in  the 
Canadian  Shield,  Bull.  Seism,  Soc.  Am.  , 53,  1,  167-210,  January. 

Evernden,  J.  F.  , 1953,  Direction  of  Approach  of  Rayleigh  Waves  and  Related 
Problems  (I),  Bull.  Seism,  Soc.  Am.  , 27,  393. 

Metz,  W.  D.  and  A.  L.  Hammond,  1974,  Geodynamics  Report:  Exploiting  the 
Earth  Sciences  Revolution,  Science,  183,  735,  February  22. 

McEvilly,  T.  V.  , 1964,  Central  U.  S.  Crust  - Upper  Mantle  Structure  from 
Love  and  Rayleigh  Wave  Velocity  Inversion,  Bull.  Seism.  Soc.  Am.  , 

54,  6,  1997  - 2015,  December. 

Rodriguez,  R.  G.  , 1969,  Atlas  of  Asia  and  Eastern  Europe  to  Support  Detec- 
tion of  Underground  Nuclear  Testi  ig;  Volume  V:  Crust  and  Mantle 
Conditions,  Department  of  the  Interior,  U.  S.  Geological  Survey, 

F ebruary. 

Santo,  T.  and  Y.  Sato,  1966,  World-Wide  Survey  of  the  Regional  Character- 
istics of  Group  Velocity  Dispersion  of  Rayleigh  Waves,  Bull,  of  the 
Earthquake  Inst.  , 44,  939-964. 


/ 


VII- 1 


Tryggvason,  E.  1965,  Dissipation  of  Rayle-gh  Wave  Energy,  J.  Geophys. 
Res.,  70,  6,  1449-  1455,  March  15. 

Turnbull,  L.  S.  , D.  F.  D.  Sun,  and  J.  S.  Shaub,  1973,  Determination  of 

Seismic  Source  Parameters  from  Frequency  Dependent  Rayleigh  and 
Love  Wave  Radiation  Patterns,  Semi-Annual  Technical  Report  No.  1, 
Texas  Instruments  Inc.,  Dallas,  Texas,  November  15. 


VlI-2 


